Dissipative relativistic fluid dynamics: a new way to derive the equations of motion 

from kinetic theory 
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We re-derive the equations of motion of dissipative relativistic fluid dynamics from kinetic theory. 
In contrast to the derivation of Israel and Stewart, which considered the second moment of the 
Boltzmann equation to obtain equations of motion for the dissipative currents, we directly use the 
latter's definition. Although the equations of motion obtained via the two approaches are formally 
identical, the coefficients are different. We show that, for the one-dimensional scaling expansion, 
our method is in better agreement with the solution obtained from the Boltzmann equation. 
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Dissipative relativistic fluid dynamics is an effective 
theory to describe the long-wavelength, low-frequency 
dynamics of various systems, with important applications 
in relativistic heavy- ion collisions and astrophysics 1]. 
However, the derivation of dissipative relativistic fluid 
dynamics from the underlying microscopic theory is not 
yet completely established. 

It can be rigorously shown that, in the non-relativistic, 
classical, dilute-gas limit, the Boltzmann equation be- 
comes equivalent to the microscopic BBGKY equations 
0. In the relativistic case and/or for quantum fluids, a 
rigorous proof does not exist, but it is commonly assumed 
that the Boltzmann equation is a reliable approximation 
to the underlying microscopic dynamics, in particular for 
dilute gases. Then, fluid dynamics can be systematically 
derived by introducing an appropriate coarse-graining 
scheme. 

The Chapman-Enskog expansion [sj is the most com- 
mon method to extract the fluid-dynamical equations of 
motion from the Boltzmann equation. However, this 
method is not suitable for relativistic systems, since it 
will inevitably lead to relativistic Navier-Stokes theory 
which displays intrinsic problems such as acausality and 
instabilities [4|-|y|. 

Israel and Stewart (IS) derived relativistic fluid- 
dynamical equations that do not exhibit this problem, 
by extending the method proposed by Grad for non- 
relativistic systems Q- In Grad's original work, the 
single-particle distribution function is expanded around 
its local equilibrium value in terms of a complete set 
of Hermite polynomials However, the generalization 
of Grad's approach to relativistic systems is non-trivial, 
since it is not easy to find a suitable set of orthogonal 
polynomials which could replace the Hermite polynomi- 
als 0, [ifl] • Thus Israel and Stewart introduced another 
approximation, the so-called 14-moment approximation 
[ll|, where the distribution function is expanded as a 
Taylor series in momentum around its local equilibrium 
value. The expansion is truncated at second order in mo- 
mentum and only 14 coefficients remain to describe the 



distribution function. 

It is important to note that the derivation of Israel and 
Stewart contains one additional approximation besides 
the 14 moments ansatz: they used the second moment of 
the Boltzmann equation to extract the equations of mo- 
tion for the dissipative currents and, hence, to determine 
the transport coefficients [§-[ll|. However, this choice to 
extract the equations of motion is ambiguous, because 
any moment of the Boltzmann equation will lead to a 
closed set of equations, once the 14-moment approxima- 
tion is applied. The transport coefficients appearing in 
the final equations depend on the choice of the moment. 

Thus, the choice of the moment is quite an impor- 
tant issue. In fact, it was confirmed that, at least for 
some cases, the IS equations are not in good agreement 
with the numerical solution of the Boltzmann equation 
[l3, [3- Also, the transport coefficients obtained by Is- 
rael and Stewart do not coincide with quantum-field the- 
oretical calculations [l3| ■ These inconsistencies may arise 
because of an inappropriate choice of the moment equa- 
tion. 

Then, which moment should be used to derive fiuid dy- 
namics? Remember that we are interested in the equa- 
tions of motion for the dissipative currents. Further- 
more, the dissipative currents are well-defined in terms of 
the single-particle distribution function from the kinetic 
point of view. Therefore, we can calculate the equations 
of motion for all the dissipative currents directly from 
their definitions without referring to an arbitrary mo- 
ment of the Boltzmann equation. The purpose of this 
letter is to derive new fluid-dynamical equations follow- 
ing this idea. We shall show that the form of these equa- 
tions is the same as in IS theory, but the values of the 
coefficients are different. For the one-dimensional scal- 
ing expansion, we demonstrate that the new equations 
agree better with a numerical solution of the Boltzmann 
equation than the IS equations. 

We start from the relativistic Boltzmann equation 



K^d^fK^C[f], 



(1) 
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where K^^ = (E'k, k) with i?k = vk^ + m^ with m being 
the particle mass. In the cohision term we consider only 
elastic two-to-two collisions, 

C[I] ^ \j dK'dPdP' Wkk,^pp, 

X (fpfp'fKfK' - fKfK'fpfp') . (2) 



Here, dK = gd^K/ [(27r)'^i?k] is the Lorentz-invariant 
measure, with g being the degeneracy factor, and 
WfCKf^ppi is the transition rate of the collision. We used 
the notation fx = fix^\K^') and /a' = 1 - af{x'',Kt'), 
where a = 1 (a = — 1) for fermions (bosons) and a = 
for a Boltzmann gas. 

The conserved particle current N'^ and the energy- 
momentum tensor T'"^ are expressed in terms of the 
single-particle distribution function as 



ATA' ^ (ifM) ^ 



(3) 
(4) 



where (...)= / dK {. . .) f, 



K- 



We introduce the fluid four- velocity u^^ as an eigenvec- 
tor of the energy-momentum tensor, T^'^Up, = eu'^ , where 
the eigenvalue e is the energy density 15[ . Then, we can 
decompose the four-momentum as 



K^' = {u- K) li^ + if 



(5) 



Here, we defined the scalar product of two four-vectors 
A^, as Afj^Bf^ = A ■ B and we introduced the pro- 
jection operator A''" = g^''' - u^'u'' and A^''^ = A^'^A^ 
for an arbitrary four- vector A^. The metric tensor is 



diag(-|-. 



Using this decomposition, A^'' and T^'^ can be written 
in the form, 



^f^^ (p + n) 



(6) 



where the particle density n, the particle diffusion current 
nf^, the energy density e, the shear stress tensor tt'^", and 
the sum of thermodynamic pressure, Pq, and bulk viscous 
pressure, H, are defined by 



= /k^^'k''^) , Po + n = - ^ 



where A^f""^ = A^'''°'^Aa0, with A'^""^ = 
(A'^^A'^" + a^'^a/'m _ |A^''A"'5)/2. We de- 
fine the local equilibrium distribution function as 
foK = (exp {/3o u ■ K — ao) -I- a) ^, where /3o and ao are 
the inverse temperature and the ratio of the chemical 
potential to temperature, respectively. These are defined 
by the matching conditions 



no = (u • K)a ,e = eo = {{u- Kf)^ 



where (. . .)o = / dK (. . .) fox- 

The separation between thermodynamic pressure and 
bulk viscous pressure is then achieved by 

^0 = {A'^-'K^K,)„ , n = (A^'-X^if.), , (9) 

with (...), ^(. ..)-(. 

So far, there is no difference between the calculation 
of Israel and Stewart and ours. The difference emerges 
in the derivation of the equations of motion for the dissi- 
pative currents. Israel and Stewart obtained these equa- 
tions from the second moment of the Boltzmann equation 
Mil] 

(K'^K'^K^) = J dK K^K^ C [/] . (10) 

Then, the equations of motion for H, n^, and tt^'^ 
are obtained by the projections u^uxd^l^K'^K'^K^'^^ 
Alu^d^{Kf'K''K^), and A"f 9^ (X^/sT^X^), respec- 
tively, together with the 14-moment approximation for 
the single-particle distribution function (see below). 
These equations determine the time evolution of H, n^, 
and TT^" through their comoving derivatives, H, g^^^ = 
A^'''q^, and 7t<^''> ee A'"'"^7rQ/3, respectively, where A = 
u ■ OA is the comoving derivative. 

However, we can calculate these comoving derivatives 
also directly from Eq. ([7]): 

ri = -im^ / dK6f , (11) 

(12) 

TT^f-' = / dA A A () r • (13) 



n<''> = J dKK'^f'Uf , 
^ J dKK^^'K^Hf 



Then, using the Boltzmann equation ([T]) in the form 

5f = -/o 



-i-if .V/+^C[/] , (14) 
u ■ K u ■ K 



where V = A^'^d^,, we obtain the exact equations 
ri = -C - /3n^ ~ CnnHfl + Cn.vr^V^, 

-Cn„a.n + ^V,((«.if)-iX<''))^ 

9 
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+ —{{u-K)-'K'^K'')^V^u,, (15) 

+ho (Hu'^ - V^n) + hoA^^'dxT:^ 
-Aii\7^(^(u- Ky^K^^^K^"^ 

- (^{u- K)-^K^''>K°'K^'j^V^Ui3, (16) 
3 



-A^^Vp ({u ■ K)-^K''^^K^''^K^p'> 



(8) 



(^{u ■ K)-^K'^>' K"^ K"Kl^^ ^"^I ^up, (17) 



3 



where ho = no/{eo + Po), and we introduced the vorticity 
WAp = I (VaUp - VpMA), the shear tensor ctap = V^A^p) 
and the expansion scalar 6 = V^u^. Above, we used the 
following notation for the collision terms, 

C j dK{u-K)-'C[f] , 

^ J dK{u- Ky^ K''''^C[f] , 

C^"" = y" dK{u-Ky^ K^^'K^^Cif] . (18) 

However, because the remaining terms in angular 
brackets cannot be entirely expressed in terms of the 
macroscopic variables (O, Eqs. ([T5l) . (fT6|). and (fT7|) are 
not closed. In order to obtain a closed set of equations, we 
use the 14-moment approximation for the single-particle 
distribution function introduced by Israel and Stewart 

Ik = foK + foKfoK (Ann + A„n„/f" + X^iTo^pK^K^) , 

(19) 

and insert this into Eqs. ([T5|) . (IT51) . and ([T7)) to compute 
the terms in angular brackets. This system of equations is 
now closed, since the approximation (jl9p solely involves 
the quantities of Eq. The coefficients An, A„ and 

A^ are well-known functions of u- K, ap, and /Sg, see e.g. 
Refs. [i-ini for details. 

We finally obtain the equations of dissipative relativis- 
tic fluid dynamics. 



ri = -— - /^n^* - ^nn9 • n - Tn„n • it - (5nnII6' 
m 



-Xnnn ■ Van + AnTrvr^'^cTui/ 



(20) 



-A„„nV^; + A„nnV^ao - A„,^''''V,ao , (21) 



(22) 



The derived equations contain 25 transport coefficients, 
of which we only show the following three coefficients 
explicitly. 



/3n 



\ - c?) (eo + ^o) - ^ (eo - SPq) 
■-^(("•^) >o, 



(23) 



Pn = ^{l), + ^{{u ■ K)-\ - ^ho, (24) 

OPO OPO PO 

= ^Po + ^(eo-3Po)-^((^i•i^)~')o,(25) 

where the velocity of sound (squared) is cj. = 
{dPo/deo)sg/no where sq is the entropy density. The other 



coefficients will be reported in Ref. ^Td\. While the form 
of the derived equations (piU) , (PT|) , and (1^ are the same 
as those obtained in previous calculations [13, the 
transport coefficients are different. That is, the deriva- 
tion of the equations of dissipative relativistic fluid dy- 
namics from the Boltzmann equation is ambiguous and 
depends on the method applied. We remark that, in the 
non-relativistic (low-temperature) limit, the set of trans- 
port coefficients as computed with the method of Israel 
and Stewart and ours converge to the same values. 

Now we would like to quantify the difference between 
the IS equations and ours at hand of a simple exam- 
ple. We consider a massless Boltzmann gas equation 
of state and the one-dimensional Bjorken scaling expan- 
sion, where the velocity is given by = i(t,0, 0,z), 
and the fluid-dynamical variables are only a function 
of the proper time t. Then, 11 and vanish. The 
shear stress tensor tt^'' has only diagonal components 
and can be characterized by a function it as tt^"^ = 
diag(0,V2,7r/2,~^).^ 

The equation for tt is given by 



dir TT ^4 , TT 

-^ + — =/3.:^-A-. 
dr Ttt or r 



(26) 



In the massless limit, our transport coefficients simplify. 



(3. 



5 ' 



■fxPoA, A 
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"63" 



(27) 

where a is the total cross section [Igj]. Here we assumed 
that cr is independent of energy and momentum as is done 
in Refs. [l^, ll^l- As mentioned above, the form of Eq. 
pT]) is identical to that of IS theory, but the transport 
coefficients assume different values. In IS theory, these 
coefficients are given by 



2Po 
3 ' 



= -aFo/3o, A = 2 



(28) 



Equation (|26|) couples to the equation of the pressure 
which is given by 



dr 3r 3t 



(29) 



In Fig. [U we show the time dependence of the 
anisotropy of the effective pressure, which is defined by 



Pl{t) _ Po{t)-tt{t) 



Pt{t) Po(t) -f ^(t)/2' 



(30) 



We used T = 500 MeV and tt = as initial condition. 
The solid and dashed lines represent our result and the 
result of IS theory, respectively. The circles correspond 
to the numerical solution of the Boltzmann equation 13]. 
This calculation is performed with values for the cross 
section such that the shear viscosity 77 = Pt^Tt^ to en- 
tropy density s ratio is constant. Since all the results 
must be compared by fixing a common cross section, 
■q = 4/(3cr/3o) and 77/5 = 6/(5(t/3o) (the shear viscosity 
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FIG. 1: Time evolution of the ratio Pl/Pt for our equations 
(solid line), the IS equations (dotted line), and the numerical 
solution of the Boltzmann equation (circles) . 

of the IS theory) have different values and are related by 

One can see that IS theory (the complete IS equations 
[l8| l always overestimates the anisotropy obtained by the 
numerical solution of the Boltzmann equation, even for 



very low viscosities {rjis/s = 0.05). On the other hand, 
our equations clearly show a better agreement. Visible 
deviations are only observed for the case oirj/s = 3.33 at 
late times. This result indicates that our fluid-dynamical 
approach is better adapted than IS theory to capture the 
microphysics contained in the Boltzmann equation. 

In summary, wc have proposed a new method for de- 
riving the fluid-dynamical equations from kinetic theory. 
In our approach, the equations for the dissipative cur- 
rents are obtained directly from the definitions of these 
currents. This method is different from the traditional IS 
approach [11], where the equations are extracted from the 
second moment of the Boltzmann equation. Our method 
can successfully reproduce the numerical solution of the 
Boltzmann equation for the simple one-dimensional scal- 
ing expansion. It is also important to mention that the 
transport coefficients of our kinetic calculation are con- 
sistent with those calculated from quantum field theory 
with the method proposed in Ref. [lj|. 
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